LPT-Orsay-01/58 
BNL-NT-01/9 



Heavy Quark Colorimetry of QCD Matter 

Yu.L. Dokshitzer a and D.E. Kharzeev 0,5 



°) LPTHE, Universite Pierre et Marie Curie (Paris VI), 
4, Place Jussieu, 75252 Paris, 
and 

LPT, Universite Paris Sud, 91405 Orsay, France 

b ) Physics Department, 
Brookhaven National Laboratory, 
Upton, New York 11973-5000, USAQ 

Abstract 

We consider propagation of heavy quarks in QCD matter. Because 
of large quark mass, the radiative quark energy loss appears to be qual- 
itatively different from that of light quarks at all energies of practical 
importance. Finite quark mass effects lead to an in-medium enhance- 
ment of the heavy-to-light D/ir ratio at moderately large (5-10 GeV) 
transverse momenta. For hot QCD matter a large enhancement is 
expected, whose magnitude and shape are exponentially sensitive to 
the density of colour charges in the medium. 
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1 Introduction 



The energy loss of heavy quarks in QCD matter has excited considerable 
interest. In particular, collisional energy loss in quark-gluon plasma has 
been evaluated by the methods of finite-temperature QCD (see for a 
comprehensive review). However, the energy loss of fast partons in medium 
is dominated by radiation of gluons and, to the best of our knowledge, a 
consistent treatment of radiative energy loss of heavy quarks in QCD matter 
is still lacking. Practical importance of this issue has become clear from a 
number of studies which showed a strong dependence of charmed hadron 
and lepton spectra on the assumed heavy quark energy loss. 

In this letter we follow the recent analysis of the quenching of inclusive 
particle spectra in QCD medium (7j to compare the yields of hadrons from 
light and heavy quarks produced in hard interactions in heavy ion collisions. 
We find that the specific pattern of gluon bremsstrahlung off heavy quarks 
(the dead cone effect) results in a much smaller heavy quark quenching, 
leading to a relative enhancement of heavy particle production in heavy ion 
collisions. 



2 Medium induced gluon radiation — a brief 
summary 

We first recall the basic features of gluon radiation caused by propagation of 
a fast parton (quark) through QCD medium. 

As was pointed out in j4|, the accompanying radiation is determined by 
multiple rescattering of the radiated gluon in the medium. The gluon, during 
its formation time 

10 . . 

tf orm — ~^2~ ' (-U 

accumulates a typical transverse momentum 

k 2 ± ~ /x 2 (2) 
A 

with A the mean free path and fi 2 the characteristic momentum transfer 
squared in a single scattering. This is the random walk pattern with an 
average number of scatterings given by the ratio tf orm /X. 
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Combining (|2]) and (|I|) we obtain 



describing the number of scattering centres which participate, coherently, in 
the emission of the gluon with a given energy u. For sufficiently large gluon 
energies, u > p 2 X, when the coherent length exceeds the mean free path, 
N co h > 1. In this situation the standard Bethe- Heftier energy spectrum per 
unit length describing independent emission of gluons at each centre gets 
suppressed: 



dW ( dW \ BH 1 _ a s C R fpJx 



duidz \duidz J N co h ttlu X V w 




(4) 



Here Cr is the "colour charge" of the parton projectile (Cr — Ci 
4/3 for the quark case we are interested in). 



2N C 

In @ we have substituted the characteristic ratio /i 2 /A by the so-called 
gluon transport coefficient || 

^ P \% q2dq ^ (5) 

which is proportional to the density p of the scattering centres in the medium 
and describes the typical momentum transfer in the gluon scattering off these 
centres. 

The transport coefficient for cold nuclear matter was expressed in as 

^I^H^g 2 )], (6) 

with p ~ 0.16 fm -3 the average nuclear density and the gluon density 

in a nucleon. The latter should be evaluated at the characteristic scale Q 2 ~ 
qL, which makes strictly speaking, an equation for q. Given a slow 
logarithmic Q 2 dependence of the gluon density, this equation can be solved 
iteratively. Choosing for the sake of estimate L ~ 5 fm results in Q ~ 
0.5 GeV. Taking at this scale a s ~ 0.5 and [xG(x)j ~ 1 (at x < 0.1), yields 

g cold ~ 0.01 GeV 3 ~ 8 p. (7) 
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This value is an agreement with the result of the analysis of the gluon p± 
broadening from the experimental data on J/ if) transverse momentum distri- 
butions |J, which in the present notation yielded 

q= (9.4 ±0.7) p. (8) 

An estimate || for a hot medium based on perturbative treatment of gluon 
scattering in quark-gluon plasma with T ~ 250 MeV resulted in the value 
of the gluon transport coefficient of about factor twenty larger than ([/]): 

g hot ~0.2 GeV 3 ~ 20q cold . (9) 

Multiplying (H) by the length L of the medium traversed,^ we arrive at the 
following expression for the inclusive energy distribution of gluons radiated 
by a quark: 

dW a s C F [uoj _ 2 / in x 

' — , uj < uji = qL . (10) 



duo ixu) V uj 

A more accurate treatment shows that the semi-quantitative estimate (|j) is 
valid for uj < cji/2, while at very large gluon energies uj ^> uj\ the gluon yield 
is very small and decreases as dl/duo oc co> -3 (for details, see ||). 

The fact that the medium induced radiation vanishes for uj > uj\ has a 
simple physical explanation, as according to (^) the formation time of such 
gluons starts to exceed the length of the medium: 



" form 



A • 



/ UJ 

/i 2 A 




Another important feature of medium induced radiation is the relation 
between the transverse momentum and the energy of the emitted gluon. 
Indeed, from ([!]) and @) (see also @) we derive 

(11) 

This means that the angular distribution of gluons with a given energy uj is 
concentrated at a characteristic energy- (and medium-) dependent emission 
angle 

/„ / * \ V4 

(12) 




2 For the sake of simplicity we assume here that the medium is static and uniform. 
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3 Radiation off heavy quarks: dead cone 



Gluon bremsstrahlung off a heavy quark differs from the case of a massless 
parton (produced in a process with the same hardness scale) in one respect: 
gluon radiation is suppressed at angles smaller than the ratio of the quark 
mass M to its energy E. Indeed, the distribution of soft gluons radiated by 
a heavy quark is given by 

a s C F duj k\dk\ n- M n Q \ 

where the strong coupling constant a s should be evaluated at the scale deter- 
mined by the denominator of ([13]). Equating, in the small-angle approxima- 
tion, k± with uj9 we conclude that the formula (|13"D differs from the standard 
bremsstrahlung spectrum 

,„ a s Cpduidk 2 , a s Cpdojd9 2 ,„ lS 

dP ~ = ^^—-^ 14 

7T UJ kj_ 7T UJ 9 2 

by the factor 

dP m = dP • (l + |) (15) 

This effect is known as the "dead cone" phenomenon. Suppression of 
small-angle radiation has a number of interesting implications, such as per- 
turbative calculability of (and non-perturbative A/M corrections to) heavy 
quark fragmentation functions |^0|,P|, multiplicity and energy spectra of 
light particles accompanying hard production of a heavy quark ]12[ . 

In the present context we should compare the angular distribution of 
gluons induced by the quark propagation in the medium with the size of 
the dead cone. To this end, for the sake of a semi-quantitative estimate, we 
substitute the characteristic angle (0) into the dead cone suppression factor 
( |T5D and combine it with the radiation spectrum ([|) to arrive at 



where 



dW a s C F RUT 1 

I{u) = uj— — = J — — xo m9 , (16) 

v ' duo 7T V uo (1 + (£uj) 3 / 2 ) 2 v ' 



-1/3 (M\^ 



r 1 " ( i ) • (17) 
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The suppression factor fllCf ) is displayed in Fig. [T] for two quark energies 
(x = u/p±). 

To see whether the finite quark mass essentially affects the medium in- 
duced gluon yield, we need to estimate the product £u for the maximal gluon 
energy u ~ u)\ to which the original distribution @ extends: 

ft* = <r" 3 (f ) V3 • ^ = (3s) 4/3 , £„ Q s m^. (is) 

This shows that the quark mass becomes irrelevant when the quark energy 
exceeds the characteristic value .Ehq which depends on the size of the medium 
and on its "scattering power" embodied into the value of the transport coef- 
ficient. 

Which regime is realized in the experiments on heavy quark production 
in nuclear collisions? Taking M = 1.5 GeV for charm quarks and using the 
values (0) and (|9|) we estimate 

/ / L \ 3/2 

Ehq = v^coid £ 3/2 M ~ 20 GeV (— ) , (19) 
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E m = V^ot L 3 / 2 M ~ 92 GeV ( —J , (20) 

for the cold and hot matter, respectively. We observe that for the transverse 
momentum (energy) distributions of heavy mesons the inequality E -C -Ehq 
always holds in practice, especially for the hot medium. We thus conclude 
that the pattern of medium induced gluon radiation appears to be qual- 
itatively different for heavy and light quarks in the kinematical region of 
practical interest. 



4 Quenching 



The issue of in-medium quenching of inclusive particle spectra was recently 
addressed in 0. The p± spectrum is given by the convolution of the trans- 
verse momentum distribution in an elementary hadron-hadron collision, eval- 
uated at a shifted value p± + e, with the distribution D(e) in the energy e 
lost by the quark to the medium-induced gluon radiation: 



da 



med 



dp] 



de D( 



da v 



dp] 



■(p± + e) = 



da v 



dp] 



■(p±) • Q(Pj 



(21) 



with Q(p±) the medium dependent quenching factor. The two facts, namely 
that in the essential region e < pi and that the vacuum cross section is a 
steeply falling function, allow one to simplify the calculation of the quenching 
factor Q by adopting the exponential approximation for the e- integral in (|21|): 



Q(px)- JdeD(e) expj— -A, C 



d 



In 



da 1 



dp] 



■(Pj 



p± J d\np± 
This integral results in the Mellin moment of the quark distribution 

C 



(22) 



D(v) = exp 



—v 



did N(lu) e 



P± 



(23) 



where N(u) is the integrated gluon multiplicity defined according to (see [0 
for details) 
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Here we present an approximate evaluation based on a simplified energy 
spectrum (|16|). Integration of flTBD leads to the following expression for the 
gluon multiplicity 



N{u) 



2a s C F 



(25) 



where 



K i x ) = ~r 



2tt 



x 



2 . 
+ - In 
9 



3^ 3(1 
1 — Jx + x 



z 3 / 2 ) 3v/3 



arctan 



(1 + v^) 5 



(26) 



The function K(x) has essential support in the region £ou = x -C 1, that 
is where the gluon yield is not suppressed by the quark mass (dead cone) 
effect (cf. fll6f)). To find the quenching factor we need to evaluate its Mellin 
transform according to (p3"|) , with v given in (^). The essential energies uj 
in the multiplicity function in ( 23|) are then restricted to < 1/za It is easy 
to check that, due to the bias effect (large value of C ~ 10), this restriction 
is stronger than the quark mass energy cutoff: 



Q 



£x4/3 



M) 



C 



M 4 



1/3 



> 2 



(27) 



for p± > 5 GeV, even for cold matter. Therefore, in this practically interest- 
ing domain, we can approximate 



N(u) 



2a s C F 



and get 



DM 



71 



cxp 



1 



9VS 



2a s C F 



7T 



qL 



} [ 1 + O (£uj 
9^3 ) 



TTV — 



(28) 



(29) 



The use of 



Qh(p 



furnishes our final result: 

2a s C F 



± 



exp 



r LCh , 16a s C F ( q M 2 



1/3' 



(30) 
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The first term in the exponent in (^) represents the quenching of the trans- 
verse momentum spectrum which is universal for the light and heavy quarks, 



Ql(p±) - exp 



2a s C F L LLl 
i- ( pi 



(modulo the difference of the C parameters determined by the p± distri- 
butions in the vacuum). The second term is specific for heavy quarks. It 
has a positive sign, which means that the suppression of the heavy hadron 
p± distributions is always smaller than that for the light hadrons. This is a 
straightforward consequence of the fact that the heavy quark mass suppresses 
gluon radiation. At very high transverse momenta, both terms vanish - this 
is in accord with the QCD factorization theorem, stating that the effects 
of the medium should disappear as p± — > oo. How fast this regime is ap- 
proached depends, however, on the properties of the medium encoded in the 
value of the transport coefficient q and in the medium size L. 

Constructing the ratio of the quenching functions, we estimate the heavy- 
to-light enhancement factor as 



Qh(p±) 
Ql(p±) 



exp 



16a s Cp 




(31) 



This simple expression provides a reasonably good approximation to more 
accurate quantitative results presented below. 



5 Quantitative estimates 

To perform a quantitative calculation of the transverse momentum distribu- 
tions, our considerations have to be combined with a realistic treatment of 
nuclear geometry, including its centrality dependence, gluon multiple inter- 
actions in the initial state (the "Cronin effect"), and a model for the time 
evolution of QCD matter in the final state. Such an analysis is beyond the 
scope of the present paper. Nevertheless, we provide some semi-quantitative 
illustrations of the expected consequences of our results, based on a simplified 
model of a static uniform medium and a fixed quark path length L. 

In what follows we use the full inclusive soft gluon emission spectrum off 
a massless quark from ||, taking into account the first non-soft correction 
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0(x): 

Hx) 



1 — x) F h In f cosh 2 



7T 



sm 




a; 



(32) 



which, for the heavy quark case, is supplied with the dead-cone suppression 



factor ([151) . To illustrate the "dead cone" effect, in Fig. |2] a comparison of 
the inclusive one-gluon emission spectra off light and charm quarks is shown, 
for hot medium with L = 5 fm. 



1.5- \ 




x 



Figure 2: Comparison of energy distributions y / x/(x) of gluons radiated 
off charm (solid line) and light (dashed line) quarks in hot matter with 
q = 0.2 GeV 3 (p ± = 10 GeV, L = him). 

As we have discussed above, the quenching of p± distributions of heavy 
hadrons caused by QCD matter is much weaker than for pions. This is 
because the gluon cloud around the heavy quark is "truncated" by the large 
quark mass, in a way which depends on the properties of the medium. This 
interesting effect can be illustrated by the transverse momentum dependence 
of the ratio of hadrons originating from the fragmentation of heavy and light 
quarks (for example, D/tt ratio) in heavy ion collisions. 

To estimate this ratio we need to know the behaviour of the vacuum 
spectra. For light hadrons (which we assume to originate from light quarks^) 

3 we remark that the quenching of hadrons produced in the fragmentation of primary 
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the parameterization of the p± distribution which describes the first RHIC 
hadroproduction data for p± up to ~ 6 GeV |H| : 



— = A ( 1 Y (33) 
dp\ \Po + P± J 

with po = 1.71 GeV and m = 12.42. 

The p± distributions of D mesons produced in hadron collisions were 
experimentally found [|14]] to be well described by the following simple pa- 
rameterization: 

d ( 1 \ ^ 

dp\ ^ \bm 2 c +p\J ' ^ ^ 

with 6 = 1.4 ±0.3, n = 10.0 ±1.2, and m c = 1.5 GeV. This parameterization 
also provides a good fit to the theoretical calculations \TE\ in perturbative 
QCD. Using this parameterization in fl22|) gives 



v = — T , a = bm r . (35) 

a 2 + p z ± 

for the v parameter in (|23|). 

Fig. |3] shows the ratio of quenching factors ( |30"D for heavy and light quarks 
in cold nuclear matter (L = 5fm), which is relevant for large-p_i_ particle 
production in pA collisions. A small value of the transport coefficient, g co id — 
0.01 GeV 3 , translates into a ~ 15% enhancement of the heavy-to-light ratio 
(about 1% for L = 2fm). 

According to (|31|), we expect a larger quenching ratio for a hot medium. 
Indeed, as Fig. |4] demonstrates, the D /it ratio should become significantly 
enhanced as compared to pp collisions: assuming a fixed length L = 5 fm of 
the hot medium traversed by the quarks, we find a factor of ~ 2 enhancement 
at p ± ~ 5 -T- 10 GeV. 

We only present the ratio of quenching factors for the following reason. 
As noticed in 0, in the region of transverse momenta we are considering, 
the absolute magnitude of quenching turns out to be extremely sensitive to 
gluon radiation in the few-hundred- MeV energy range and therefore cannot 
be quantitatively predicted without detailed understanding of the spectral 
properties of the medium. 

gluons contains an additional factor 9/4 ~ 2 in the exponent 
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H - L Ratio 

1.25r 




Figure 3: The ratio of quenching factors Qh(p±)/Ql(p±) for charm and light 
quarks in cold nuclear matter (q = 0.01 GeV 3 , L = 5fm). 

The heavy-to-light ratio, however, proves to be much less sensitive to the 
infrared region, since gluon radiation off heavy and light quarks is universal 
in the x — > limit, see fll6l). To illustrate this point, in Fig. (| we show the 
ratio of quenching factors calculated with the energy restriction upon gluon 
energies, u > 500 MeV. We see that in the 5 — 10 GeV range of p±, this 
modifies the ratio by 20 — 30 %fj 

The predicted ratio should not be taken at its face value: the enhancement 
factor originates mainly from a very strong quenching of light quark jets, and 
in reality such a dramatic suppression will be washed out by light-particle 
production at the periphery of the collision zone. In a central AA collision, 
the number of collisions in this peripheral zone, depending on the strength 
of absorption, scales as ~ A 1 / 3 A 2 ^ 3 , whereas the total number of collisions 

4 The quenching factors themselves change (increase) by an order of magnitude when 
the radiation of gluons with energies smaller than 500 MeV is vetoed. 
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which contribute to the production of D mesons scales as ~ A 4 / 3 . Thus, 
even with account of the collision geometry, we may still expect an order of 
~ A 2 / 3 -T- A enhancement of the D/ti ratio, which for a Au — Au collision 
may translate in a strong effect. 

Moreover, given a large nucleus and moderately large p±, both light and 
(especially) heavy quarks will be tempted to turn into hadrons prior to leav- 
ing the interaction arena. In the pion case, this is likely to add to quenching 



via pion absorption in the medium [ 16| , thus enhancing further the D/ir 
ratio. 

Clearly, detailed calculations have to be performed before a defendable 
number for the magnitude of the D/ir enhancement can be presented. Nev- 
ertheless, the D/tt ratio appears to be extremely sensitive to the density of 
colour charges in QCD matter, and we eagerly await the results of experi- 
mental studies of this quantity in relativistic heavy ion collisions. 
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Figure 4: The ratio of quenching factors Qh(p±)/Ql(p±) f° r charm 
and light quarks in hot matter with q = 0.2 GeV 3 (L = 5fm upper 
panel; L = 2fm lower panel). Solid lines correspond to unrestricted 
gluon radiation, while the dashed lines are based on the calculation 
with the cut on gluon energies u > 0.5 GeV. 
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